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Research  Abstract 


Objectives: 

The  primary  objective  of  this  research  is  to  implement  and  test  improved  mixed  layer 
parameterizations  in  a  layered  numerical  model  of  the  Indian  Ocean  circulation.  This 
will  allow  inclusion  of  thermodynamic  forcing  as  well  as  wind  stress  forcing  of  the 
model  circulation  in  order  to  understand  the  variability  and  complicated  mixed-layer 
physics  observed  in  the  real  Indian  Ocean  and  to  develop  an  optimum  formulation  for 
mixed  layer  physics  in  layered  ocean  models.  Although  these  efforts  are  focused  on  the 
Indian  Ocean,  die  formulation  developed  here  can  be  extended  to  other  areas  of  the 
world  ocean. 

Approach: 

The  ocean  circulation  model  is  that  described  in  Luther  and  O'Brien  (1985)  as  extended 
by  Jensen  (1991)  and  by  J.  Capella  (personal  communication).  The  present  version  of 
the  model  has  4  layers,  with  the  horizontal  pressure  gradient  assumed  to  vanish  in  the 
lowest  layer  (the  reduced  gravity  approximation).  The  model  domain  covers  the  Indian 
Ocean  basin  from  30  degrees  S  to  26  degrees  N  and  from  35  degrees  E  to  120  degrees  E 
at  a  resolution  of  1/12  degrees  in  latitude  and  longitude.  Thermodynamic  forcing 
through  a  simplified  Kraus-Tumer  mixed  layer  parameterization  as  in  McCreary  et  al. 
(1993)  has  been  added  to  this  version  of  the  model.  Entrainment  and  detrainment 
processes  are  included  that  allow  water  to  move  between  layers.  Detrainment  simulates 
the  subduction  of  surface-layer  water  into  the  deeper  ocean.  The  entrainment  and 
detrainment  processes  are  dynamically  similar  to  the  vertical  mixing  of  momentum  and 
heat  that  is  present  in  continuously  stratified  models.  This  improved  model  of  the 
Indian  Ocean  is  being  evaluated  against  available  observations  and  in  comparison  with 
a  parallel  model  formulation  being  undertaken  at  NRL  by  Kindle  and  Gallacher. 

Tasks  completed: 

A  20-year  integration,  driven  by  climatological  monthly  mean  winds  from  the  FSU 
analysis,  has  been  used  to  "spin-up"  the  present  version  of  the  model  (with  wind  stress 
forcing  only).  Using  the  end  of  this  20-year  integration  as  an  initial  condition,  a  model 
integration  has  been  completed  using  the  monthly  mean  real-time  winds  processed  at 
Florida  State  University  for  the  period  January  1977  through  February  1994.  This 
integration  is  updated  through  the  16th  of  the  most  recently  completed  month  by  the 
5th  working  day  of  each  month  to  provide  a  near-real-time  "quick-look"  estimate  of 
conditions  in  the  Indian  Ocean.  The  results  from  this  model  integration  are  available  to 


the  research  community  via  anonymous  ftp  from  kelvin.marine.usf.edu.  These  data  are 
being  used  in  planning  the  upcoming  ONR  Accelerated  Research  Initiative  in  the 
Arabian  Sea  and  WOCE  and  JGOFS  work  in  the  Indian  Ocean.  From  December  22, 

1992,  to  March  9, 1993  there  were  258  remote  entries  into  the  anonymous  ftp  directory 
from  numerous  US  and  international  internet  sites  (we  stopped  keeping  track  of  ftp 
entries  after  that  period). 

Interannual  variability  in  the  model  circulation  fields  has  been  investigated  for  the 
period  1977-1992.  A  composite  climatological  seasonal  cycle  has  been  removed  from 
the  monthly  mean  fields  of  velocity,  layer  thickness,  and  sea  surface  elevation  to  create 
monthly  anomaly  fields.  A  12-month  running  mean  has  been  applied  to  the  monthly 
anomaly  fields  to  isolate  the  interannual  variability.  A  complex  empirical  orthogonal 
function  (CEOF)  analysis  is  being  performed  on  the  monthly  anomaly  fields  to  identify 
the  dominant  modes  of  variability. 

Implementation  of  the  mixed-layer  formulation  of  McCreary  et  al.  (1993)  into  the 
circulation  model  code  has  been  completed.  Luther  and  a  Ph.D.  student,  Zaihua  Ji,  are 
now  testing  this  version  of  the  circulation  model  with  the  assistance  of  J.  P.  McCreary. 
This  model  code  has  been  given  to  John  Kindle  at  NRL-SSC  for  evaluation.  The 
accompanying  figures  show  the  model  upper  layer  thickness,  velocity,  and  mbced  layer 
temperature  for  May  15  and  August  15  after  12  years  of  spin-up  using  the  climatological 
atmospheric  forcing  fields  from  Rao  et  al.  (1989,/.  Geophys.  Res.,  94, 10801-10815).  After 
testing  is  completed,  parallel  runs  of  our  model  and  the  NRL  model  will  be  run,  using 
identical  forcing  fields.  Results  of  these  model  runs  will  be  compared  against  each  other 
and  against  observations.  - 

Under  separate  funding  from  NSF  and  NASA,  a  biogeochemical  model  has  been 
developed  and  is  being  coupled  to  the  physical  model.  The  improved  mixed-layer 
formulation  in  the  physical  model  provides  information  on  mixed  layer  depth,  nutrient 
entrainment  and  mixing  intensity  required  by  the  biogeochemical  model.  The  coupled 
physical-biogeochemical  model  will  be  applied  to  analysis  and  interpretation  of  the 
observations  from  the  ONR-JGOFS  Arabian  Sea  Process  Study. 

Results: 


Significant  interannual  variability  is  found  in  the  model  fields.  The  dominant  large- 
scale  interannual  variability  occurs  in  the  eastern  tropical  Indian  Ocean  and  has  similar 
characteristics  to  the  Pacific  ENSO  phenomena  in  that  it  is  instigated  by  eastward- 
propagating  equatorial  Kelvin  waves.  Associated  with  this  mode  of  variability  is  a 
large-scale  redistribution  of  mass  in  the  upper  layer  of  the  ocean  north  and  south  of  10 
degrees  south  latitude.  The  periods  September  1979  to  December  1981  and  April  1983 
to  June  1985  appear  as  warm  anomalies  in  the  eastern  tropical  Indian  Ocean,  with  the 
thermocline  as  much  as  15  m  deeper  than  usual.  The  periods  Februaty  1987  to  August 
1987  and  January  1988  to  December  1990  appear  as  cold  anomalies  with  a  shallower 
thermocline  in  the  eastern  tropical  Indian  Ocean.  Gary  Meyers  and  R.  R.  Rao  (personal 
communication)  see  similar  variability  in  XBT  observations  from  the  eastern  Indian 


Ocean. 


Interannual  variability  in  the  western  basin  is  of  smaller  spatial  scale  and  is  associated 
with  year-to-year  variations  in  the  Somali  Current-EACC  system  and  the  atmospheric 
Findlater  Jet.  Year  to  year  variations  in  upwelling  and  primary  production  in  the 
Arabian  Sea  during  the  period  of  1979  through  1982  as  presented  in  Brock  and  McClain 
(1991)  can  be  understood  in  terms  of  changes  in  the  strength  of  the  wind  stress  curl  field 
associated  with  the  Findlater  Jet.  Changes  in  the  wind  stress  curl  field  drive  changes  in 
Ekman  pumping,  leading  to  changes  in  the  depth  of  the  mixed-layer,  the  thermocUne, 
and  the  nutricline.  These  changes  are  clearly  seen  in  the  variability  in  thermocline 
depth  in  the  interannual  model  integration  described  above,  leading  to  a  better 
understanding  of  the  relationship  between  Ekman  pumping  and  mixed-layer  evolution 
as  modelled  by  Bauer  et  al.  (1991). 

Accomplishments: 

The  most  important  accomplishment  to  date  from  this  work  is  the  implementation  of 
the  model  in  near-real-time  and  the  distribution  of  the  results  to  the  research 
community  for  use  in  plaiming  future  ONR7TGDFS,  and  WOCE  observations  in  the 
Indian  Ocean.  By  running  the  model  with  observed  winds  from  the  FSU  analysis,  we 
have  provided  a  greater  understanding  of  the  variability  in  the  Indian  Ocean 
circulation. 

Another  significant  accomplishment  is  the  implementation  of  the  mixed-layer 
parameterization  in  the  Indian  Ocean  circulation  model.  The  model  code  has  been 
delivered  to  John  Kindle  at  NRL-SSC.  Continued  development  and  testing  of  the  model 
will  be  done  in  collaboration  with  Kindle  and  McCreary. 
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Figure  legend: 

Mixed-layer  temperature  (TM),  and  upper  layer  (thermocline)  depth  and  upper  layer 
velocity  for  May  15  and  August  15  after  12  years  of  spin-up  forced  by  the  climatological 
monthly  mean  atmospheric  fields  from  Rao  et  al.  (1989,/.  Geophys.  Res.,  94, 10801-10815). 
Atmospheric  fields  consist  of  wind  stress,  wind  speed,  air  temperature,  and  humidity. 
The  model  computes  surface  fluxes  using  these  fields  and  its  own  mixed  layer 
temperature  from  bulk  aerodynamic  fomulae.  Color  represents  temperature  or 
thermocline  depth  and  arrows  represent  velocity  (see  scale  bars  and  arrows).  The 
model  reproduces  the  warming  of  the  Arabian  Sea  prior  to  the  southwest  monsoon 
onset  and  the  cooling  after  the  onset.  The  cold  upwelling  wedges  observed  in  the  two- 
gyre  system  of  the  Somali  Current  are  simulated  by  the  model. 
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